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Quantum Field Theory is plagued by divergences in the attempt to calculate physical quantities. 
Standard techniques of regularization and renormalization are used to keep under control such a 
problem. In this paper we would like to use a different scheme based on Modified Dispersion Rela- 
tions (MDR) to remove infinities appearing in one loop approximation in contrast to what happens 
in conventional approaches. In particular, we apply the MDR regularization to the computation of 
the entropy of a Schwarzschild black hole from one side and the Zero Point Energy (ZPE) of the 
graviton from the other side. The graviton ZPE is connected to the cosmological constant by means 
of of the Wheeler-DeWitt equation. 
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I. INTRODUCTION 



q-( The appearance of a trans-Planckian physics in Black Hole thermodynamics has led many authors to consider that 
some deep change in particle physics should come into play. One realization of these ideas is represented by the 
modification of the Heisenberg uncertainty relations, better known as Generalized Uncertainty Principle (GUP)[l|-lj|- 
This principle is based on the following inequality 



A 2 

AxAp>h+-^(A P ) 2 , (1) 

©: h 

£f) , where h is the Planck constant and X p is the Planck length. Of course, the above inequality affects the Liouville 
\& • measure which becomes 

d 3 xd 3 p 

(2nhf(l + \pif {2) 

When A = 0, the formula reduces to the ordinary counting of quantum states. If Eq.@ is used for computing 
the entropy of a black hole from a Quantum Field Theory point of view, the usual UV divergence at the horizon 
can be removed[l|, Indeed, without introducing Eq. (|14l) one is forced to use traditional methods for removing 
divergences: for example renormalizing the Newton constant or using Pauli-Villars regularization 0. It is clear 
that the distortion of the Liouville measure plays a key role in regularizing divergent integrals. Non-commutative 
geometry provides another powerful method to have such a distortion. As shown in 043 , one finds [lOj 

al 3 xd 3 k , d 3 xd 3 k ( 9 j2 \ 

an — =- =>• an = 5- exp — k . ( o ) 

(2n) 3 (2tt) 3 V 4 y 

This deformation corresponds to an effective cut off on the background geometry. The UV cut off is triggered only by 
higher momenta modes > 1/vO which propagate over the background geometry. When 9 — 0, the formula reduces to 
the ordinary counting of quantum states. An application of non-commutative geometry to the computation of black 
hole entropy shows that the usual horizon divergence disappears [llj. In connection with these ideas, in recent years, 
there has been a proposal on how the fundamental aspects of special relativity can be modified at very high energies. 
This modification has been termed Doubly Special Relativity (DSR)[12]. In DSR, the Planck mass is regarded as an 
observer independent energy scale. One of its effects is that the usual dispersion relation of a massive particle of mass 
m is modified into the following expression 

E 2 g\ (E/Ep) - p 2 g\ (E/E P ) = m 2 , (4) 
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where g\ (E/Ep) and gi (E/Ep) are two unknown functions which have the following property 

lim g l (E/Ep) = l and lim g 2 (E/E P ) = l. (5) 

E/Ep^O E/Ep^O 

Thus, the usual dispersion relation is recovered at low energies. Eqs. (|4l 5| are a representation of 11 Modified Dispersion 
Relations" (MDRs). The common motivation in using them is in that they can be used as a phenomenological approach 
to investigate physics at the Planck scale, where General Relativity is no longer reliable. Moreover, we expect the 
functions g\ (E/Ep) and gi (E/Ep) modify the UV behavior of quantum fields in the same way as GUP and Non- 
commutative geometry do, respectively. Note that GUP and MDR modifications are strictly connected fllij. Since 
the form of gi (E/Ep) and 172 (E/Ep) is unknown and they have to obey the property ([5]), we have a large amount 
of arbitrariness in fixing the dependence on E/Ep, even if some specific choices have been proposed by G. Amelino- 
Camelia et al.[l4| in the context of black hole thermodynamics. MDRs play a relevant role also when the background 
is curved. Following the analysis of Magueijo and Smolin[l5| one can define the following "rainbow metric" 

N 2 (r)dt 2 dr 2 r 2 ,~ . 2 . , l7 . 

91 (E/Ep) (i_m\gl(E/E P ) 9l (E/Ep) 



which is a solution of the distorted Einstein's Field equations 

(E) = 8ttG (E) T m „ (E) + g^A (E) . (7) 

G (E) is an energy dependent Newton's constant, defined so that G (0) is the physical Newton's constant. Similarly we 
have an energy dependent cosmological constant A (E) . The function b (r) will be referred to as the "shape function" . 
The shape function may be thought of as specifying the shape of the spatial slices. If the equation b(r w ) = r w is 
satisfied for some values of r, then we say that the points r w are horizons for the metric ([5]). In this paper we will 
discuss the use of MDRs on black hole entropy calculation [l6j and on the estimation of the cosmological constant 
computed with the help of a revisited Wheeler-DeWitt equation (WDW) [13] • Units in which h = c = k = 1 are used 
throughout the paper. 

II. BLACK HOLE ENTROPY WITH MDRS 

To start with, we fix the ideas on a real massless scalar field described by the action 

I = -\ ( d*x^g-[g^d^d v ct>] (8) 



2 . 

in the background geometry of Eq. © with N (r) described by 



A 2 (r)=exp(-2A(r))(l-^l, (!)) 



where A (r) is known as the "redshift function" that describes how far the total gravitational redshift deviates from 
that implied by the shape function. Without loss of generality we can fix the value of A (r) at infinity such that 
A (00) = 0. The Euler-Lagrange equations are 

1 =9 M (V=g<Td v ) = 0. (10) 



In order to use the WKB approximation, if <p has a separable form, we can define an r-dependent radial wave number 
k(r,l,E) 



k 2 (r,l,E)= ' 



1 _ 



e, V (2A(r)) E2h2 ^ l « + 



1 Hr) 



(11) 



with 



(12) 
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The number of modes with frequency less than E is given approximately by 



(21 + 1) / y/k 2 (r,l,E)drcU, 



(13) 



where it is understood that the integration with respect to r and / is taken over those values which satisfy r w < r < R 
and k 2 (r, I, E) > 0. Thus, from Eq.(Jll|) we get 



d~g(E) 
dE 



dis(l,E) 
dE 



(21 + l)dl 



2 d (\ 



ix dE \ 3 



E 6 h 6 (E) 



(14) 



Defining f3 as the inverse temperature measured at infinity, the free energy is given by 

1 



& Jo 



ln(l 



e ., E) d9jE) dE= r 
' dE 



(15) 



We divide the integration range into two intervals: in [r^ri] we define i^and in [ri,+oo) with r\ > r w we define 
Fr. Assuming that A (r) < oo, Vr e [r w , +oo), Fr is dominated by large volume effects for large R. This term will 
give the contribution to the entropy of a homogeneous quantum gas in flat space at a uniform temperature T . We fix 
our attention on 



Fr 



2 1 ' x 

7T /? Jo 



2 exp(3^l(r)) 



(16) 



In proximity of : 



1 — = [l-o {r w )\ 



(17) 



and the radial part of F rw becomes divergent. This ultraviolet divergence has been cured by 't Hooft, who introduced 
a "brick wall tq" proportional to Zp[18(. Nevertheless, since spacetime is modified by a "rainbow metric" , it is quite 
natural that even the "brick wall" is affected by this distortion. To see such an effect, we perform the radial integration 
in F r to obtain 



drr 



2 exp(3A(r)) 



2 exp(3^(r)) 4 exp(3A(r,„)) 1 
drr 7T — r„ 



r(E/E P ) (l-^) 2 W {l-V{r w )fr(E/E P y 



(18) 



where we have assumed that, in proximity of the throat the brick wall is no longer a constant but it becomes a function 
of E/Ep. Plugging Eq.(TT8]) into Eq.® we find 



F r 



2rj exp (3A (r w )) 
ZpK(i-b> [r w )f Jo 



oo . 
31,3'- " 



E 6 h s (E/Ep) 



dE 



In (1- exp 



r (E/Ep) 



dE, 



(19) 



where we have integrated by parts with the condition that h (E/Ep) be chosen in such a way to allow the convergence 
when E/Ep — > oo. Without loss of generality we write 



r (E/Ep) = r w a (E/E P ) . 



(20) 



with 



a (E/Ep) 0, E/Ep^Q. 



(21) 



In this way the horizon divergence is still present but it is translated in terms of a function of E/Ep. Plugging Eq. ([2H|) 
into Eq. (fl9|) . we obtain 



F r ... = - 



C r 



3/3r 



w Jo 



E^(E,E P) - 



ln(l - exp 



a (E/Ep) 



dE 
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where we have defined 

c = 2r^, exp (3 A (r w )) 

It is clear that g\ (E/Ep) and 52 (E/Ep) must be chosen in such a way to compensate the vanishing of a (E/Ep), oth- 
erwise the horizon divergence (brick wall) cannot be eliminated. For example, one good candidate for the convergence 
is 

h(E/E P )=exp(-^-y (23) 
A good candidate, but not exhaustive for a (E/Ep) is 

a (E/Ep) = (-^-) , a > 0. (24) 



E P J 



In the limit where j3Ep ^> 1, the total energy U is 



d{fiF r J _ 2 exp(2A{r w )) 2E 2 P 2 exp (2 A (r w )) E% 

OP w l-b'(r w ) 9/3 2 k w 11 w l-b'(r w ) 9/3 [ ' 



and the entropy S is 



l2 dF Tw _ 2 exp(2A(r w )) 1E$ _ A Tw E 2 P exp (2 A (r w )) 2 



P S/3 w 9? Kw 4 1 - 6' (r w ) Qtt' 1 J 

where we have used the expression for the surface gravity in the low energy limit 

k w = ^— exp (-A (r w )) [1 - b' (r w )] (27) 

and where we have integrated over E. To recover the area law, we have to impose that 

exp(2A(r„)) _ 9n 
1 - V (r w ) = ~2 

and 



(28) 



?- r -s- (29) 

This corresponds to a changing of the time variable with respect to the Schwarzschild time. The total energy becomes 

U = rl^f, (30) 
which in terms of the Schwarzschild radius r w — 2MG and inverse Hawking temperature (3 — 8ttMG becomes 

U = 4M 2 G 2 -^- = — . (31) 

16A/G 4 v ' 

Note the discrepancy of a factor of 3/2 with the 't Hooft result. 



III. GRAVITY'S RAINBOW AND THE WDW EQUATION 



The WDW equation was originally introduced by Bryce DeWitt as an attempt to quantize General Relativity in a 
Hamiltonian formulation. It is described bvfl9| 



W& = 



(2k) G l3 kiK 



v^(3#_2A) 



* = 



(32) 
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and it represents the quantum version of the classical constraint which guarantees the invariance under time 
reparametrization. Gy/rf is the super-metric, ir^ is the super-momentum, 3 R is the scalar curvature in three di- 
mensions and A is the cosmological constant, while k = 8irG with G the Newton's constant. In this way, the WDW 
equation is written in its most general form. The main reason to use such an equation to discuss renormalization 
problems is related to the possibility of formally re-writing the WDW equation as an expectation value computation. 
Rather than reproduce the formalism, we shall refer the reader to Refs.[20| for details, when necessary. However, for 
self-completeness and self-consistency, we present here a brief outline of the formalism used 1 . Multiplying Eq. (|32[) by 
'J* [gij] and functionally integrating over the three spatial metric g^ we find 



1 / V [ 9ij ] [ 9ij ] J v r/ :i .r.\v ( l' fry] 1 (* k ' / '''' A 



*» A 



V /7V (/ ,, vl;- , A , v|, !h ; V (vj, 1 Of) 

In Eq. (|33[) we have also integrated over the hypersurface £ and we have defined 



V= [ d 3 x^g. 



(33) 

K 



(34) 



V is the volume of the hypersurface £ and 

A E = (2k) G ljkl ^ir kl - V? 3 i?/ (2k) . (35) 

In this form, Eq. ((33|) can be used to compute ZPE provided that A/k be considered as an eigenvalue of As. In 
particular, Eq.([3"3"f represents the Sturm-Liouville problem associated with the cosmological constant. To solve Eq. (|3"3")l 
is a quite impossible task. Therefore, we are oriented to use a variational approach with trial wave functionals. The 
related boundary conditions are dictated by the choice of the trial wave functionals which, in our case are of the 
Gaussian type. Different types of wave functionals correspond to different boundary conditions. The choice of a 
Gaussian wave functional is justified by the fact that ZPE should be described by a good candidate of the "vacuum 
state" . To fix ideas, we choose the line element (j6)) as background metric with g\ (E/Ep) = g^ (E/Ep) — 1, namely 
MDRs do not distort the metric. Then we consider a perturbation of the metric tensor of the form g^ = g^ + hij, 
where gij is the background metric and hij is a quantum fluctuation around the background. Thus Eq. (|33j) can be 
expanded in terms of hij. Since the kinetic part of As is quadratic in the momenta, we only need to expand the 
three-scalar curvature J d 3 x^/g 3 R up to the quadratic order. To proceed with the computation, we need an orthogonal 
decomposition on the tangent space of 3- metric deformations [22j| : 



The operator L maps £j into symmetric tracefree tensors 

(LOij 



V,< y + - - gij (V • • 



(36) 



(37) 



hfj is the traceless-transverse component of the perturbation (TT), namely g^hfj — 0, V 4 /i^ = and h is the 
trace of hij. It is immediate to recognize that the trace element a = h — 2 (V • £) is gauge invariant. The same 
decomposition can be done also on the momentum it ij and induces the following transformation on the functional 
measure Vhij — > Vh^V^iDaJi, where J\ is the Jacobian related to the gauge vector variable £j. The only physical 
information is encoded 

(2) 



1 



A- 



V 



After having functionally integrated, we find 



(\&-L|*-L) 



1 



ijkl 



(2k) K- 11 - (x, x) m + — (A L ) K 1 x) 



(2k) 



i hi 



(38) 



(39) 



1 See also Ref. |2l( for an application of the method to a / (R) theory. 
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where 



= (A^ x ) - AR<thi 3 + 3 ifc£. 



is the modified Lichnerowicz operator and A^is the Lichnerowicz operator defined by 



(A L h) %1 = Ah^ 



2R lk]l h kl + R lk h) + R jk hf A = -V a V Q . 



Qijki re p resen ts the inverse DeWitt metric and all indices run from one to three. Note that the term —ARfh^j 
disappears in four dimensions. The propagator K (x, x) iakl can be represented as 



(40) 

(41) 
- •"'/.•///, 



iakl 



E 



2A(r) 



(42) 



where h < £ >± (~af ) are the eigenfunctions of A L . t denotes a complete set of indices and A (t) are a set of variational 
parameters to be determined by the minimization of Eq. ([39|) . The expectation value of A^ is easily obtained by 
inserting the form of the propagator into Eq. (|39p and minimizing with respect to the variational function A (r) . Thus 
the total one loop energy density for TT tensors becomes 



A 
8^G 



(43) 



The above expression makes sense only for uif (r) > 0, where uii are the eigenvalues of A L . For a background of the 
form of Eq.©, if we define the reduced fields /j (x) — Fi (x) /r, we find that the Lichnerowicz operator (^A,h A 
be reduced to 



d 2 1(1 + 1) 2 , ; 
— + 2 + m 2 (r) 



dx 



f i (x)=uhf i (x) i = l,2 



(44) 



where we have used the Regge- Wheeler representation|23j and m 2 (r) and m| (r) are two r-dependent effective masses 

nl(r) = £(l- b -&) + £*b>(r)-£sb(r) 

^W = ^(l-^) + ^'(r) + ^6(r) 



(r = r (x)) . 



(45) 



In Eq. (|44|) , we have defined 



dx = ± 



dr 



1 



b(r) 



(46) 



Like in the entropy calculation of section (JTTJ) , we use the W.K.B. method and we define two r-dependent radial wave 
numbers 



k 2 (r, I, u^) = u 2 nl - l JL±H _ m 2 ( r ) i = 1,2 



For every degree of freedom of the graviton we apply Eq. (fl"4|) and we find that Eq. (|43]) becomes 

2 r+oc 

T I 

7T 



A 
8^G 



1=1 



1) 



dg{u)i) 

uii — duJi. 

duji 



(47) 



(48) 



This is the one loop graviton contribution to the induced cosmological constant. The explicit evaluation of Eq. (|4"8")l 

gives 



8ttG 



Pi + P2 = 



I z r+oc . 



(49) 
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where we have included an additional Att coming from the angular integration. p\ and p 2 are divergent and traditionally 
the use of the zeta function regularization keeps the divergences under control. To this purpose we reconsider Eq. (|32[) 
in presence of Gravity's Rainbow and we find 2 



*> A 



v W) " (50) 



where 



As = (2k) 3 G ljk in J TT --7=—, , , (51) 

5-2 {E/Ep) (2k) g 2 {E/E P ) 

The symbol "~" indicates the quantity computed in absence of rainbow's functions g\ (E/Ep) and g 2 (E/Ep). Of 
course, Eqs. ([5fJl) and (I5T1) reduce to the ordinary Eqs. (l32l 33[) and (f3"51 when E/Epi —> 0. By repeating the procedure 
leading to Eq. (|39|) . we find that the total one loop energy density becomes 



A i z r+°° a \ f f]2 \ 3 



8irG 



where E* is the value which annihilates the argument of the root. To further proceed, we choose a form of g\ (E/Ep) 
and g 2 (E/Ep) suggested by a Noncommutative geometry analvsisflOj. If we fix 

3 i (E/Ep) = (l + (3^j exp(-a|J) and g 2 (E/E P ) = 1, (53) 

with a > and j3 £ R, Eq. (|52p can be easily integrated. However, it is more useful to give the asymptotic expansion 



for large and small x, where x = \Jrn\ 2 (r) /Ep. The asymptotic expansion for large x is 

_A (2/3a 3 / 2 + 0Fa 2 ) x 8/3a 5 ^ + 3y7a 3 3 16/3a 7 / 2 + 5y^a 4 

8^G ~ 4a 7 /2 ieoTV^ + 128 a 15 / 2 x 3 

while for small x, one gets 

A _ 4a 5 / 2 + 3v^a 2 
8ttG ~ 4a 9 / 2 

If we set 



O (-T- 4 ) , (54) 



O (x 6 ) . (55) 



0=-^, (56) 

then the linear divergent term of the asymptotic expansion (|54j) disappears. This means that A/8irG vanishes for 
large x. On the other hand for small x we get 

A 37r — 8 _ , o N 

+ O (x 3 ) , (57) 



8ttG 8a 2 



where we have used the result of expansion (|55[) . It is possible to show that with choice (|56| , the induced cosmological 
constant is always positive. Note that when /3 = 0, the pure "Gaussian" choice can not give a positive induced 
cosmological constant [l^. 
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